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We study the middle cerebral artery blood flow velocity (MCAfv) in humans using transcranial 
Doppler ultrasonography (TCD). Scaling properties of time series of the axial flow velocity averaged 
over a cardiac beat interval may be characterized by two exponents. The short time scaling exponent 
(STSE) determines the statistical properties of fluctuations of blood flow velocities in short-time 
intervals while the Hurst exponent describes the long-term fractal properties. In many migraineurs 
the value of the STSE is significantly reduced and may approach that of the Hurst exponent. This 
change in dynamical properties reflects the significant loss of short-term adaptability and the overall 
hyperexcitability of the underlying cerebral blood flow control system. We call this effect fractal 
rigidity. 
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Migraine headaches have been the banc of humanity 
for centuries, afflicting such notables as Ceasar, Pascal, 
Kant, Beethoven, Chopin, and Napoleon. However, their 
aetiology and pathomechanism have not to date been sat- 
isfactorily elucidated. Herein we demonstrate that the 
scaling properties of the time series associated with cere- 
bral blood flow (CBF) significantly differ between those 
of normal healthy individuals and migraineurs. The re- 
sults of our data analysis show that the complex patho- 
physiology of migraine frequently leads to hyperexcitabil- 
ity of the cerebral flow control system. 

Physiological signals, such as CBF time series, are typ- 
ically generated by complex self-regulatory systems that 
process inputs with a broad range of accessible values. 
Even though this type of time scries may fluctuate in an 
irregular and complex manner, they frequently exhibit 
self-affine or fractal properties which can be character- 
ized by a single global parameter - the fractal dimension 
D or equivalently the Hurst exponent H (D = 2 — 
The salient property of mathematical random fractal 
process is the existence of long-range correlations for 
H 7^ 1/2. The studies of the cardiac beat-to-beat vari- 
ability have shown the existence of strong long-range 
correlations in healthy subjects and demonstrated the 
breakdown of correlations in disease (see also and 
references therein). A similar pattern was observed in 
fluctuations in the stride interval in human gait. The 
strength of correlations was significantly reduced both 
by aging and a neurodegenerative disease. This effect is 



frequently referred to as the loss of complexity 
Complexity decreases with the convergence of the Hurst 
exponent on H = 1/2. 

A healthy human brain is perfused by blood flowing 
laminarly through the cerebral vessels providing brain 
tissue with substrates such as oxygen and glucose. It 
turns out that CBF is relatively stable with typical values 
between 45 and 65 ml/lOOg of brain tissue per second, 
despite variations in systemic pressure as large as 100 
Torr. This phenomenon is known as cerebral autoreg- 
ulation and has been thoroughly documented not only 
in humans but also in animals 0. Autoregulation, 
which is mainly associated with changes in cerebrovas- 
cular resistance (CVR) of small precapillary brain arter- 
ies, is only one of at least four major mechanisms that 
regulate CBF. A considerable body of evidence suggests 
that CBF is influenced by local cerebral metabolic ac- 
tivity. As metabolic activity increases so does flow and 
vice versa. The actual coupling mechanism underlying 
this metabolic regulation is unknown but most likely it in- 
volves certain vasoactive compounds (which affect the di- 
ameter of cerebral vessels) such as adenosine, potassium, 
prostaglandines which are locally produced in response to 
metabolic activity. External chemical regulation is pre- 
dominantly associated with the strong influence of CO2 
on cerebral vessels. An increase in carbon dioxide arterial 
content leads to marked dilation of vessels (vasolidation) 
which in turn boosts CBF while a decrease produces mild 
vasoconstriction and slows down CBF. The impact of the 
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sympathetic nervous system on CBF is often ignored but 
intense sympathetic activity results in vasoconstriction. 
This type of neurogenic regulation can also indirectly af- 
fect cerebral flow via its influence on autoregulation. 

The complex cerebral flow regulation mechanisms are 
supposed to be influenced or even to be fundamentally 
altered in many pathological states. However, despite 
the significant advances in brain diagnostic imaging tech- 
niques many functional aspects of CBF regulation are not 
fully understood 01- F° r example, migraine - a preva- 
lent, hemicranial (asymmetric) headache is among the 
least understood diseases. In migraine without aura at- 
tacks may involve nausea, vomiting, sensitivity to light, 
sound, or movement. These associated symptoms distin- 
guish migraine from ordinary tension-type headaches. In 
about 15% of migraineurs headache is preceded by one 
or more focal neurological symptoms, collectively known 
as the aura. In migraine with aura the associated symp- 
toms may include transient visual disturbances, march- 
ing unilateral paresthesias and numbness or weakness in 
an extremity or the face, language disturbances, and ver- 
tigo 0,^3- According to the leading hypothesis, mi- 
graine results from a dysfunction of brain-stem or dien- 
cephalic nuclei that normally modulate sensory input and 
exert neural influence on cranial vessels, see, for example, 
[l2l Hii ] and references therein. Thus, the fundamental 
question arises as to whether migraine can significantly 
influence cerebral hemodynamics. Some experimental 
data reveal clear interhemispheric blood flow asymme- 
try in some parts of the brain of migraineurs even during 
headache-free intervals [Til H5I . 

Transcranial Doppler ultrasonography enables high- 
resolution measurement of MCA flow velocity. Even 
though this technique does not allow us to directly de- 
termine CBF values, it may help to elucidate the na- 
ture and role of vascular abnormalities associated with 
migraine. Some previous studies have shown significant 
changes in cerebrovascular reactivity in migraine patients 
[l7j . In this work we look for the signature of the mi- 
graine pathology in the scaling properties of the human 
MCAfv time series. 

The dynamical aspects of the cerebral blood flow reg- 
ulation were recognized by Zhang et al. [isj . Keuncn et 
al. 0,H3| applied the attractor reconstruction technique 
along with the Grassberger-Procaccia algorithm and the 
concept of surrogate data to look for the manifestations of 
the nonlinear dynamics in continuous waveforms of TCD 
signals. Rossitti and Stephensen used the relative 
dispersion of the MCAfv velocity time series to reveal 
its fractal nature. West et al. |22| extended this line of 
research by taking into account the more general prop- 
erties of fractal time series. Both studies |2l|,|22j showed 
that the beat-to-beat variability in the flow velocity has 
a long-time memory and is persistent with the average 
value of the Hurst exponent H = 0.85 ± 0.04, a value 
consistent with that found earlier for interbeat interval 



time series of the human heart. Finally, West et al ob- 
served that cerebral blood flow is multifractal in nature 

H 

We measured MCAfv using the Multidop T DWL Elek- 
tronische Systeme ultrasonograph. The 2-MHz Doppler 
probes were placed over the temporal windows and fixed 
at a constant angle and position. The measurements were 
taken continuously for approximately two hours in the 
subjects at supine rest. The study comprised 15 healthy 
individuals and 33 migraineurs (14 had migraine with 
aura and the others had migraine without aura). Mi- 
graine was diagnosed according to the guidelines of the 
Headache Classification Committee of the International 
Headache Society [24| ■ An example of a typical measured 
MCAfv time series is shown in Fig. \I\for the first thou- 
sand of the recorded beats of the subject's heart. The 
total time series has over eight thousand data points for 
a two hour data record. 

Successive increments of mathematical fractal random 
processes are independent of the time step. They are cor- 
related with the coefficient of correlation p which is deter- 
mined by the formula 2 2H = 2 + 2p. Thus for H 1/2 
there exist long-range correlations, that is, p ^ 0. It 
turns out the Hurst exponent also determines the scaling 
properties of the fractal time series. If y{t) is a frac- 
tal process with Hurst exponent H, then y c — y(ct)/c H 
is another fractal process with the same statistics. The 
variance of self-affine time series is proportional to At 2H 
where At is the time interval between measurements. A 
number of algorithms which are commonly used to cal- 
culate the Hurst exponent are based on this property. 

Herein we employ the detrended fluctuation analysis 
(DFA) introduced into the study of biomedical time series 
by Peng et al. |25| . Let {vi}^ =1 be the experimental time 
series of the middle cerebral artery blood flow velocity 
(MCAfv) v. First, the time series is aggregated: y(k) = 
Yli=i( v i ~ v), k = 1, .., N, where v is the average velocity. 
Then, for segments of the aggregated time series of length 
n the following quantity is calculated: 
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where y no is a least square line fit to the data segment 
which starts at no and ends at n + tiq. The bar in the 
above equation denotes an average over all possible start- 
ing points no of data segments of length n. Thus, for a 
given data box size n, F(n) gives the characteristic size of 
fluctuations of the aggregated and detrended time series. 
If the aggregated time-series is fractal then F{n) ~ n H , 
so one obtains the Hurst exponent from a linear least- 
square fit to F(n) on double log graph paper. However, 
West has emphasized 0, |2^ the importance of possible 
periodic modulations of quantities such as F(n). These 
modulations may be accounted for with the help of the 
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following fit function: 

F x (n) = n H exp[a + Xcos(j\nn)]. (2) 

Here again the Hurst exponent is determined by the slope 
of the fitting curve, but now the curve also has a harmonic 
modulation in the logarithm of the length n of the data 
segment. 

Fig. |3 shows the typical DFA analysis for a healthy 
subject. The circles in this figure are the calculated val- 
ues of F(n) and the solid line is the best renormalization 
group fit, cf. Eqn. J5J). The best-fit parameters are given 
in the inset. It is apparent from this graph that frac- 
tal properties of the cerebral flow, as indicated by the 
quality of the fit, are well pronounced only for time in- 
tervals larger than approximately 32 cardiac beats. Scal- 
ing properties for shorter intervals are distinctly different 
and may be characterized by a short time scaling expo- 
nent (STSE). The slope of the grey line in Fig. H yields 
the value of the STSE. The average value of the STSE 
for the control group was 1.34 ± 0.11 while the average 
Hurst exponent was H = 0.80 ± 0.10. The averaging was 
done over 36 calculations such as the one shown in Fig. 
121 We would like to point out that the values of both 
exponents are not affected by a scaling of the time series 
amplitude and do not depend on the series's mean value 

cf. ©. 

To elucidate the nature of the two scaling regions char- 
acteristic of healthy individuals let us consider the follow- 
ing one-dimensional map: 

Vi+i =Vi- bvi + — + £i, (3) 

which models the fluctuations of the blood flow velocity 
averaged over a cardiac beat. The map is reminiscent of 
the Mackey- Glass differential equation originally intro- 
duced to describe production of white blood cells [26| . 
The role of the linear term in the above equation is to 
dampen out fluctuations. The nonlinear positive feed- 
back term with delay is the source of long-range correla- 
tions. £ is a random variable which mimics the stochastic 
component of the cerebral hemodynamics. In this work 
we choose £ to be normally distributed with zero mean 
and standard deviation a. We call the difference equa- 
tion © a cerebral blood flow map (CBFM). Fig. |3]shows 
the DFA of a time series generated by the CBFM with 
the parameters b = 0.10, a = 0.05, r = 15 and a = 0.20. 
The map clearly exhibits the characteristics of the exper- 
imental data including the distinct crossover. 

It turns out that the scaling properties of the MCAfv 
time series may be profoundly influenced by the migraine 
pathophysiology. In our study, for about 40% of mi- 
graineurs the value of the STSE is significantly reduced 
and may closely approach the value of the Hurst expo- 
nent cf. Figs. 0] and 03 On the other hand, for both the 
migraineurs with aura and without the average Hurst 



exponent is the same as that of the control group. The 
CBFM provides the insight into dynamical origin of this 
effect. The value of the STSE is dependent mainly on 
the strength of the damping b. With increasing b the 
STSE decreases. However, to maintain approximately 
constant value of the Hurst exponent the strength of the 
positive-feedback term a must also increase. Fig. [H] ex- 
emplifies this behavior (the time series was iterated with 
the following parameters: b = 0.45, a = 0.35, r = 15 and 
<7 = 0.20). Thus, the reduction of the STSE observed in 
the experimental data is the result of the excessive damp- 
ening of the cerebral flow fluctuations and is the mani- 
festation of the significant loss of adaptability and over- 
all hyperexcitablity of the underlying regulation system. 
We call this novel effect fractal rigidity. We would like 
to emphasize that hyperexcitabity of the cerebral blood 
flow control system seems to be physiologically consis- 
tent with the reduced activation level of cortical neurons 
observed in some transcranial magnetic stimulation and 
evoked potential studies 0, 0, 0, and references 
therein. 
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FIG. 1: MCAfv time series for a healthy subject. 



10 



8 - 



Se- 



ra 4 
O H 

















STSE 


= 1.39 
















a= 
H= 


-0.8 
0.85 - 
0.23 




■ ^ 






X= 








r= 


0.97 



6 

Log 2 (n) 



10 



12 



FIG. 2: Detrended Fluctuation Analysis of the MCAfv time 
series of a healthy subject. The slope of the the grey line 
gives the value of the short time scaling exponent STSE. The 
parameters of the fit are shown in the right bottom corner 
of the graph. 




FIG. 3: DFA of a time series generated by the cerebral blood 
flow map (CBFM) with the following parameters: b = 0.10, 
a — 0.05, t = 15 and a — 0.20 (the time series was scaled to 
facilitate comparison with the DFA shown in Fig. |5J. 
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FIG. 4: DFA of the MCAfv time series of a migraineur with- 
out aura. The time series was recorded during a headache-free 
interval. 




FIG. 5: DFA of the MCAfv time series of a migraineur with 
aura. The measurement was performed during a headache. 




FIG. 6: DFA of a time series generated by the CBFM with 
the following parameters: b = 0.40, a = 0.35, r = 15 and 
a = 0.20. 



